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We investigate the quasinormal modes for gravitational perturbations of rotating black holes in 
four dimensional Anti-de Sitter (AdS) spacetime. The study of the quasinormal frequencies related 
to these modes is relevant to the AdS/CFT correspondence. Although results have been obtained 
for Schwarzschild and Reissner-Nordstrom AdS black holes, quasinormal frequencies of Kerr- AdS 
black holes are computed for the first time. We solve the Teukolsky equations in AdS spacetime, 
providing a second order and a Pade approximation for the angular eigenvalues associated to the 
Teukolsky angular equation. The transformation theory and the Regge-Wheeler-Zerilli equations 
for Kerr-AdS are obtained. 
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I. INTRODUCTION 

D 

The holographic principle in field theory states that, in the presence of gravity, physics in a region of space-time 
can be encoded on a lower-dimensional surface. The Maldacena AdS/CFT conjecture Q provides explicit examples 
for some of these holographic encodings. For example, IIB string theory on the background AdS$ x S 5 or M— theory 
on the background AdS± x S 7 are claimed to be mathematically equivalent (dual) to certain quantum field theories 
£S) , on the boundary at infinity of the relevant AdS space. The isometry groups of the AdS spaces become conformal 
' groups acting on the boundary quantum theories. In the case of AdS^ x S 5 , the quantum field theory is Af — 4 Super 
^vO \ Yang-Mills gauge theory with gauge group U(N) and coupling constant gyu- This particular Yang-Mills theory is 
\J ■ also a conformal field theory (CFT). 

More specifically, if one introduces a generating function for correlation functions of operators Oj in the Yang-Mills 
theory, then the AdS/CFT correspondence states that 0,1a], 

o : ^-f0<°>cM*)^ T = z string [<t> k ]. (i.i) 

In this equation, on the left-hand side, we have the generating function with operators integrated against arbitrary 
?H . density functions, <f>y , on the boundary of AdS. The boundary has dimension D = 4, 3 according to the version 
of the correspondence we consider. The function Zstringl&k] is the string partition function with string coupling g s , 

evaluated in the AdS background. The fields, satisfy the boundary conditions of approaching (f>y at infinity. 
A limit can be taken to reduce the string theory to its supergravity limit. This may be done by a combination 
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■ of taking the string tension to infinity {a 1 — > 0) and sending g s N — > oo. A second limit is obtained by making the 
approximation (g s — > 0) of the supergravity partition function to exp(— S), where the supergravity action S[(j>k] is 
the action of the classical supergavity solution that satisfies the boundary conditions mentioned above. Although 
the strongest form of the correspondence is the string theory expression, in practice it is much easier to work with 
supergravity actions, and it is for this reason that most results have been obtained in this context. On the quantum 
field theory side of the correspondence, this approximation amounts to taking both TV and gy M N large. The classical 
supergravity action should holographically capture the large- N limit of strongly coupled quantum Yang-Mills theory. 

A quite large amount of work has been produced in the last few years on the AdS/CFT conjecture and an 
extensive introduction can be found in 0,0- In the AdS/CFT correspondence, the study of AdS black holes has 
a direct interpretation in terms of the dual conformal field theory on its boundary. The duality predicts that the 
retarded CFT correlation functions are in one to one correspondence with Green's functions on Anti-de Sitter space 
with appropriate boundary conditions Wi ITol. Illf. Furthermore, as mentioned in |12|. it is assumed that a large 
static black hole in AdS space corresponds to a thermal state in the CFT on the boundary. Perturbing the black hole 
with a scalar field, for instance, is equivalent to perturbing this thermal state. The perturbed system is expected, 
at late times, to approach equilibrium exponentially with a characteristic time-scale. This time-scale is inversely 
proportional to the imaginary part of the poles of the correlators of the perturbation operator. It seems that these 
relaxation time-scales are quite complicated to calculate in the CFT, therefore their computation is conveniently 
replaced by the evaluation of the quasinormal frequencies in the AdS bulk space. 

In and , the authors went beyond the scalar perturbation treated in ^1 an d they considered electromagnetic 
and gravitational perturbations. In one of the most recent works on the subject, Berti and Kokkotas |l5j confirm 
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and extend previous results on scalar, electromagnetic and gravitational perturbations of static AdS black holes and 
analyse, for the first time, Reissner-Nordstrom Anti-de Sitter black holes and calculate their quasinormal frequencies. 
Naturally, perturbations associated with other kinds of fields than the scalars treated in [lj] will decay at different 
rates given by their quasinormal mode frequencies. 

Quasinormal modes have recently become relevant to quantum gravity for a different reason. This arises from 
Hod's proposal to describe quantum properties of black holes from their classical oscillation spectrum [Tfil |. The quest 
for a consistent quantization of black holes was initiated many years ago by Bekenstein |l7l lisj who, based on semi- 
classical arguments, conjectured that the horizon area of a (non extremal) black hole should have a discrete eigenvalue 
spectrum. The eigenvalues of this spectrum were expected to be uniformly spaced. Hod noted that the real parts of 
the asymptotic quasinormal frequencies of a static black hole of mass M, as computed by Nollert [l^, can be written 
as lur — Tln3, where T is the black hole Hawking's temperature. He then exploited Bohr's correspondence principle, 
asserting that 'transition frequencies at large quantum numbers should equal classical oscillation frequencies', to 
conjecture that variations in the black hole mass due to quantum processes should be given by AM = huiR. Finally, 
using the first law of black hole thermodynamics and the equivalence between entropy and the surface area, he deduced 
the equal spacing in the area spectrum of a static hole. 

Hod's conjecture has brought further motivations in the field of black hole physics and has provoked several attempts 
to evaluate highly damped quasinormal frequencies of Kerr black holes, ft is, in fact, straightforward to ask whether 
this conjecture is valid for more general black holes than the static ones; if asymptotic frequencies depend in these 
more general cases on black hole charge and/or angular momentum; if the presence of a cosmological constant should 
imply a modification of the proposal. At the moment these ideas of universality in Hod's proposal have not been 
confirmed 20)]. 

A brief outline of this paper is the following. In section^we review some of the main features of AdS spacetime and 
its static and rotating black hole solutions. In section lTTTI we solve the radial Teukolsky equation. Then, we concentrate 
on the angular equation and propose both a second order and a Pade approximation of the associated eigenvalues. In 
section IV we introduce transformation theory and obtain the Regge- Wheeler equation. These results are exploited 
for the explicit determination of the quasinormal frequencies and some considerations on their asymptotic spacing. 
All the numerical results are reported in section We conclude by a final discussion of the results we have obtained 
and with some ideas for future work. In the appendix we give the explicit recurrence relation for the series solution 
of the radial Teukolsky equation. 



II. ADS SPACETIME AND ITS ROTATING BLACK HOLE SOLUTIONS 

Rotating black holes in four dimensions with a cosmological constant were first constructed by Carter [2^] and 
Demianski [2|j. In the asymptotically de Sitter case, the black hole solutions have cosmological horizons. Thermal 
equilibrium is only possible when the black hole horizon has the same temperature as the de Sitter horizon [2j|. The 
asymptotically Anti-de Sitter case has no cosmological horizons and the black hole, even if it is rotating, can be in 
equilibrium with radiation out to the conformal boundary at infinity |24l l25l |2(j . 

Anti-de Sitter space and the associated stationary black hole solutions can be described by a metric with radial 
component g rr of the form 



r 2 + a 2 cos 2 1 



(2.1) 



The polynomial A r (r) vanishes at the event horizon, which we shall denote by r+. The event horizon is also a Killing 
horizon, that is it is a surface on which a Killing field £ a become null. The surface gravity k is defined in terms of £ a 
by, 



V Q (£ b 6) = -2< a . (2.2) 

For a stationary and axisymmetric black hole with timelike and axial killing fields dt and d^, the killing vector field 
£ a is of the form 

£ = d t + (2.3) 

The constant Qh is the angular velocity of the horizon. 

An explicit expression for the surface gravity of the event horizon is given by p2| , 



A '^ + ) (24) 
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The Hawking temperature of the black hole is related to the surface gravity through the famous result, 



T = 



2tt' 



(2.5) 



The version of the AdS/ CFT conjecture which concerns us here is the correspondence between M— theory on AdS$ x S 7 
and supersymmetric conformal Yang-Mills field theory on the conformal boundary at infinity. Crucial to the AdS /CFT 
conjecture is the fact that the AdS space and its related black hole solutions are conformal to a rotating version of the 
Einstein static Universe at infinity. The angular velocity of the space at infinity is denoted by fl. The combination 
(T, f2) form a set of thermodynamic variables [23.l25l.l2aj. 
For Kerr- AdS spacetime with cosmological constant A, 



A r = (a 2 + r 2 )(l - ^Ar 2 ) - 2Mr 
3 



(2.6) 



In Fig. ^ we show the conformal compactification for the covering space of Kerr- AdS* spacetime. The conformal 
boundary is the boundary of region I. 



ring singularity 
r = 




FIG. 1: The Penrose diagram of a Kerr black hole in Anti-de Sitter space-time is shown. IH and H are the inner and event 
horizon of the black hole. Regions I, bound by vertical segments (r = oo), represent the exterior space. Regions II are the 
regions between the inner (r = r_) and the event horizon of the black hole. Regions III correspond to r < r_. The bold 
vertical segments represent the ring singularity at r = in the equatorial plane. 



The black hole event horizon is the largest real root of A r = 0, 

(r 2 + a 2 )(l + a 2 r 2 )-2A:/r = 0, (2.7) 

where a 2 = —A/3 is the inverse of the Anti-de Sitter radius. From Eq. 12.71 and the definition of the surface gravity 
12.41 we derive an expression for the temperature T, 

3a 2 r1 + (1 + a 2 a 2 )r\ - a 2 , , 

47rr+(r 2 _ + a 2 ) 

The rotation of the asymptotic space is given by 0> HE US 

a(l + a 2 r 2 _) , 

tt=—n 2.9 

a 2 + rl 
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Fig. [21 shows contours of constant k in the aa and Ma parameter space. The n = solutions are extremal black 
holes and above this line in the diagram the solutions become naked singularities. Fixing the value of the angular 
parameter a and k we have unique masses and therefore unique solutions for aa > 0.1, but up to 3 different solutions 
for aa < 0.1. 
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FIG. 2: Isotherms related to Kerr-AdS black holes are plotted in the aa and Ma parameter space. 



The transformation from the metic parameters M and a to the thermodynamic variables becomes singular where 
the Jacobian 



d(T,Q) r 2 + {3a 2 r 2 + -l)-a 2 (l + a 2 r 2 + ) 



d(r+,a) 



^{r\ + a 2 ) 2 



(2.10) 



vanishes. Vanishing of the Jacobian also corresponds to the divergence of the heat capacity at constant il. This 
critical line is plotted in Fig. |3J Another important feature of the parameter space is set by the condition aft, < 1. If 
aft > 1, then the killing field £ becomes non-timelike beyond a certain radius and the black hole radiation can not 
remain in thermal equilibrium. This allows us to divide the paramter space into the four regions A, B, C, D. Region 
A is the region with no event horizons. Regions B and C both have negative specific heat. Region D is the region of 
positive specific heat and afl < 1. 



III. QUASINORMAL FREQUENCIES OF KERR ADS BLACK HOLES 



The full Kerr AdS metric can be found in [23|, for instance, and it can be expressed as 



A p 2 , r + M ■ 2 , r 2 + [i 2 , 2 



{r 2 +n 2 ) {r 2 +(i 2 ) A r A M 

where r, 9 and <fi are spherical coordinates, with /i = etcos^. The forms uj\ and u>2 are 

lo\ = dt — ax 2 sin 2 9 dtp, Ll>2 = dt — a _1 x _2 (r 2 + a 2 ) d<p 
and the metric functions are given by 

A r = ( a 2 + r 2 )(l - ^Ar 2 ) - 2Mr, 



(3.1) 



(3.2) 



(3.3) 
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FIG. 3: This figure of the aa and Ma parameter space shows various stability limits, afi = 1 and the critical line where the 
thermal capacity Cn diverges. The parameter space is divided into regions: A with only naked singularities; B where aQ, > 1 
and Cn < 0; C where aQ < 1 and Cn < and D where afl < 1 and Cn > 0. 



3 

We shall make use later of the radial tortoise coordinate, defined by 

d A r d 



A, = (a 2 -M 2 )(l + gA M 2 ), (3.4) 
X 2 = 1 + ^Aa 2 . (3.5) 



dr* (a 2 +r 2 )dr' 

Using the the definition of surface gravity given in Eq. 12.41 applied to each of the roots of A r = (two of them 
complex), we can express r* in the form 

r* = In a(r — r + ) H \na(r — r_) H In a(r — rs) H — — lna(r — fa), (3-7) 

2k + 2k- 2«3 2^3 

Note that r* e (— oo, 0), with r* — > as r — » oo as a consequence of the identity nZ 1 + nZ 1 + 1 + S3 1 = 0. (In the 
ordinary Kerr solution, r* — > 00 as r — > 00). 

The generalisation of Teukolsky's equations and their separation for asymptotically de-Sitter and Anti de-Sitter 
spacetimes, was first derived by Chambers and Moss [3JJ|. They looked for separable solutions for the Weyl scalar of 
the form 

y Q = e +luJt e ml ^R{r)Si(8) 1 (3.8) 

where Si(9) can be found in |3lj |. and I is the angular mode number characteristic of the spherical harmonics. (In 
|30j | , the time coordinate was scaled differently and the frequency was a factor of \ 2 larger than the frequency which 
we are using here.) In what follows we will drop the index I when we talk about the angular function Si. 
Following [30| . we introduce new radial and angular derivatives 

Un ~ ° r + A r + A r ' " r A r + A r ' [ ' 

K„ nA' . K„ nA'„ 

L n = d, + -^ + -^ 4 = ^--^ + -^, (3.10) 
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with K r and defined as 

K r = x 2 am+ (a 2 + r 2 )to, (3.11) 
Kn = x 2 am+ (a 2 ~ P 2 )uj. (3.12) 

The angular and radial functions S{9) and R(r), satisfy the two second order differential equations 

{L_ x A^L\ + 6lu[i - 2A/?)S = -AS (3.13) 
(£)_!A r £)} - 6iwr - 2kr 2 )R = \R (3.14) 

where A is a separatin constant. In the non-rotating case, A = (I — + 2) for 1 = 2,3,.... 

A. Radial equation 

Let us start from the radial differential eauation l3.14l which can be re- written as 

iK r A' 1\ A r / iK r A' 1\ -6iu Al , 

This equation has four regular singular points at the four roots of A r ; one at the black hole outer horizon, r + , which 
is the largest root of A r , one at the black hole inner horizon, r_, and other two complex conjugate roots, r 3 and f 3 , 
which are unphysical roots. This allows us to write A r as 

A,, = a (r — r + ) (r — r_) (r — r 3 ) (r — r 3 ) . (3.16) 

where a 2 = —A/3 is the inverse of the Anti-de Sitter radius. Next, we set x = r^ 1 , so that we can write 

A r = px - 4 , (3.17) 

with 

p = a 2 (x — x\) (x — sea) {x — x 3 ) (x — 5 3 ) . (3.18) 

Note that x\,X2,x^, x 3 are just the inverse of the roots of A r in terms of the new variable x. 
Using the same coordinate inversion we can rewrite K r as 

K r = -x- 2 q, (3.19) 

with 

q = -X 2 amx 2 + (l + a 2 x 2 )uj. (3.20) 

It is also useful, for future reference, to re-write p as 

p = a + (1 + a 2 a 2 ) x 2 - 2Mx 3 + a 2 x 4 , (3.21) 

In the new variable x the radial equation becomes 

iq p' 3\ f n iq p' 3\ 6iiv 2A 1 „ „ . . 

flx + -- - + -pa x --i + i--- --X\R = 0. (3.22) 

p p x J \ p p X J X X A J 

Rescale R as 

R = xp- x Y (3.23) 
and insert it into eauation l3.22l After dividing everything by x 2 , we obtain 

-3 x -^- A - + t\(\y^r + ^Y]-^Y-^Y-^Y = 0. (3.24) 
p x p J \x z x xp J x A p x A p xp 
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Rescale Y as 

Y = e +iu)r "Z, (3.25) 

where 

n=0 ^ ' 1J 

By writing the radial solution in this way we achieve the result of factoring out the behaviour of this function at the 
black hole event horizon r = r + (r* = — oo). Substitute Y and its derivatives into 13.241 We finally get an equation in 
the function Z that reads 

Z'' + ( 2 ^-l)z' + (^-^-^-^- 3 -(^- P -))z = 0. (3.27) 
\ p p J \xp x 2 - x A p p x \ p p J J 

The determination of the recurrence relation for the radial series solution Z can be found in appendix A of this paper. 



B. Angular equation 

The solution to the radial equation is strictly bound up with the determination of the eigenvalues A(a, a, to, lu) of 
the angular equation, which appear explicitly in the coefficients of the recurrence relation associated with the series 
solution Z of the radial differential equation. We could attempt to solve for both u> and A simultaneously, but this 
would be a time-consuming procedure. Instead, we shall attempt to determine A before we find u. To this purpose, 
we first propose a second order Taylor approximation of A for small auj. We then improve this with the new numerical 
technique, already used in the context of Kerr black holes [^3, based on a Pade approximating function for the 
separation constant. 



1. Second order approximation for A 

We proceed in analogy with Sasaki |33j|. who used a technique based on perturbation theory to evaluate the 
separation constant for Kerr black holes in asymptotically flat spacetime. Our final result agrees with Sasaki's result 
when we assume the cosmological constant to be equal zero. 

We shall only consider to = 0. Let us start by rewriting the angular equation in the form 

Kjj_ _ . /„ k u a;, 



9 M + ^ - ^ J A J 9 M - ^ + ^ ) + 6o,/i - 2A^ + A } S = 0. (3.28) 



Since m = 0, this implies for the expression 



The equation can also be written as 



with 



K fl = a 2 ujsm 2 e. (3.29) 



HS = -AS* (3.30) 



K 2 IkA 1 A 12 

H = Adl + A'd, — K' ft + A" H — — + 6ujfi — 2A M 2 . (3.31) 

First of all, it is straightforward to verify that when A = and a — we get from H the expression we expect for 
static black holes in asymptotically flat spacetime, as can be verified from |3lj . 
We then introduce the new variable u = cos#, such that /i = au. This gives 



A A' K 2 9hA' A' 2 

H= i±Ql + ^± du _ K > + A " - + qfl. -=- + GLuau - 2AaV. (3.32) 
a 1 a p AAA 
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Next, we expand H in powers of aw, 

H = J ff + awi/i+a 2 cj 2 i7 2 + C'(a 3 cj 3 ), (3.33) 



where, 



d 2 cos(9 d „ 4 
dO 1 sin 80 sin 

ff x = -4cos6> (3.35) 

■ M 4 + ;*- (1 - u y (3.36, 



We also expand and A in powers of aw as 

5= So + awS' 1 + a 2 uj 2 S 2 + <3(a 3 w 3 ), (3.37) 

A = A + au\ 1 + a 2 cu 2 X 2 + 0(a 3 w 3 ), (3.38) 

with A = (I - 1) Q + 2) . 

Now, if we plug 13.331 13.371 and 13. 381 into Eq. 13.301 we have at the second order 

(Ho + auH x + a 2 u) 2 H 2 )(So + auSx + a 2 uu 2 S 2 ) = 

— (Ao + auXi + a 2 uj 2 \2)(So + clujSi + a 2 uo 2 S 2 ). (3.39) 

Collecting the terms of the same order we obtain, 

(H + Xo)S l = (3.40) 
(H + X )S[ = -(H t + Xt)So, (3.41) 
(H + X )S l 2 = -(H t + Xi)S[ - (A 2 + H 2 )S l . (3.42) 

The functions S l are the spheroidal harmonics of spin weight s = — 2, for m = 0. 

We will be interested in the explicit expressions of the functions S l , which can be represented using Legendre 
Polynomials Pi (cos 8) , as 

S l (9) =Pi A e-Pi,e cote. (3.43) 

We should mention that at the moment the spherical harmonics are not normalised. It is convenient to normalise at 
least those we will be using. The normalised spherical harmonics are 



A2 _ 

J o — 


^s 2 


(3.44) 


S* 3 - i 
— \ 


' 240 


(3.45) 


<? 4 - * 
<->o — 1 


/ Q 4 

/so v 


(3.46) 



Since our purpose is to evaluate A up to the second order, let us start with the calculation of Xi in Eq. 

The first order correction to A is obtained by multiplying eauation l3.41l bv S l from the left hand side and integrating 
it over 9. This gives 

Xi = - S l HxSq sin8d0 = 4 / (5 ) 2 cos 9 sin 6d6, (3.47) 
Jo Jo 

but, if m = 0, this integral is always zero for any value of I as can, indeed, be verified in |33j . 
Consequently, 

Ai = 0, V/, m = 0. (3.48) 
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To obtain S[, we set 

s' = EW m = °- ( 3 - 49 ) 

We insert this expression into eauation l3.41l multiply it by Sq and integrate over 0. What we get is, 

<${*o-K)J (^o) 2 sin0d0 = -4 jf d(cos6)S l o cos 6^. (3.50) 

/ 

4 



The coefficients of the series on the left hand side, q , can be derived from 13.501 and written as 



(Z'-l)(Z' + 2)-(Z-l)(Z + 2) 



/ d(cos 6)S l Q cos 6S l (3.51) 
Jo 



Hence cj is generally non zero only for Z' = Z ± 1. 

The second order correction, A2, is obtained by multiplying equation 13.421 by Sq form the left hand side and 
integrating it over 9. We find, 

/■7T /-7T 

A 2 = -4/ d(cos6)S cos9S[ + S l {H 2 S )d(cos8). (3.52) 
Jo Jo 



But 



f d(cos0)S o cos0S< = ~ $>f ) 2 [(Z' - 1)(Z' + 2) - (Z - 1)(Z + 2)] (3.53) 
Jo 4 

and the sum is non zero only for I' = I ± 1. It follows that, 

A 2 = -( c | +1 ) 2 [Z(Z + 3)-(Z-l)(Z + 2)] 

- ( c ;- 1 ) 2 [(Z-2)(Z + l)-(Z-l)(Z + 2)]+ [ S l H 2 S l Q d(cos6). (3.54) 



If we, finally, evaluate the coefficients c\ and the integral on the right hand side of 13. 541 when Z = 2, what we obtain 
for A2 is 



2 A 10 

A2 =7^ + 21 ^ 



This gives us, for Z = 2, up to the second order in au, an expression for the separation constant A that reads 

2 10 

A = A + a 2 uj 2 \ 2 = 4+ -Ka 2 + — aV 2 . (3.56) 

A similar calculation can be repeated for any other value of Z. 

For instance, we have performed the same calculation for Z = 3. In that case, the Taylor expansion of A in the 
variable aw, up to the second order is 

2 

A = A +a 2 oj 2 \ 2 = 10 + Aa 2 + -a 2 uj 2 . (3.57) 

o 

In both cases, Z = 2 and Z = 3, we recover the flat space results for A as soon as we set A = 0, as can be checked in 

m. 



2. Pade approximation for X 

Having showed how to obtain a second order approximation for the angular eigenvalues, we would like to partially 
overcome the incompleteness of this approach by a different scheme, which is based on the construction of a Pade 
approximating function for A. The technique which we describe below has already been applied to find the quasi 
normal modes for Kerr black holes in asymptotically flat space |32| . and works well for all values of the parameter aw. 
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We will start from an analytical solution of the Teukolsky angular equation in AdS space-time. This equation is 

(L^ 1 A^L[+6cjfj,-2Ai/) S = -XS. (3.58) 

If we replace the differential operators L_i and h\ with their explicit expressions 13. 101 we can rewrite Vs. 581 as 

„ A' , / A" K' Kl A' A' 2 Qua 2 An 2 A \ 

S" + ^S'+ l-f--f-g + 2K^-£ + ^-^ + ^ )S = 0. (3.59) 

This differential equation has a structure which is easily recognisable as being very similar to the structure of the 
radial equation, with the variable \x having assumed the role of i times the radial variable r. It has four regular 
singular points and they correspond to the four roots of A M . They are /i = ±a and /i = ±a. 

As already done for the radial solution, it is possible to factorise out the leading behaviour at /j. — a by rescaling 
the angular function S as 

S = A^e+^'Z, (3.60) 
The variable which is analogous to the radial tortoise coordinate r*, is defined by 

d A, 



<i/i» (a 2 — /i 2 ) <i/i ' 

while 

(M - a) n 



(3.61) 



y a „^_^. (3.62) 
^ (-2a)" v ; 



n=0 



This allows us, through the same logic used to deal with the radial function R, to get to this simpler form for the 
angular equation 

A„Z" + {2K I1 - A') Z' + (6ujh - 2Afi 2 + A) Z = 0. (3.63) 



Next, we expand all the coefficients of 13.631 as Taylor series around /i = a and insert them into the equation itself 
beside the expression for Z and its derivatives. We obtain a four term recurrence relation of the form 

+ $ n a n + ln<in-i + 5 n a n -2 = 0. (3.64) 

Finally, we find the coefficients of 13.641 in the explicit form 

a n = (n 2 - 1) (l - a 2 a 2 ) (3.65) 
/3„ = -n (n - 1) (1 - 5a 2 a 2 ) + n (-Aaoj + 2- 10a 2 a 2 ) 

+6auj + 6a 2 a 2 + A (3.66) 
7„ = -8a 2 a 2 (n - 1) (n - 2) - 2 (n - 1) {-2auj - 12a 2 a 2 ) 

-12aw-24a 2 a 2 (3.67) 

5 n = A{n-A){n-h)a 2 a 2 . (3.68) 

This allows us to determinate the coefficients a n of the series solution by recursion. 

Equation 13.631 has critical points at the north and south poles /i = ±a. We require the solution which does not 
diverge at either pole. This solution has ao — 0, a\ — 1 and 

- ^2, (u-a) n 

Z = 2_^a n ,_ 2a)n -» as Li^-a. (3.69) 

Consequently, we have 

oo 

J2 an = 0. (3.70) 

n=0 
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FIG. 4: The real part of the Pade approximation of the angular eigenvalue A is shown to fit the data for lambda with imaginary 
ui. Note that the Pade approximation was constructed for real u. The angular parameter act — 0.24, angular mode number 
I — 2 and m = 0. We also show the improvement we achieve with our Pade function technique comparing our results to a 
second order approximation of A. 

At this point, we wrote a C+ + computer program able to evaluate values for A(aa, aw), which are identified with the 
roots of the left hand side of Eq. 13.701 numerically. Within the program, the values of A which satisfy Eq. (|3.70|1 
are found by the use of the Newton- Raphson method routine extended to the complex plane [HiJ ■ We have fit these 
numerical data for A to Pade approximating functions for real auo at a fixed value of aa, given by 



an + ai(aw) + a,2(auj) 2 + a 3 (auj) 3 
1 + bi(auj) + 62 (a^) 2 



(3.71) 



Finally, the Pade approximation and the numerical data are compared for imaginary auj. In the Pade formula, the 
numerator is one degree higher than the denominator. The reason for this choice is that the Teukolsky angular equation 
belongs to the family of differential equations for angular prolate spheroidal wave functions, and the behaviour of the 
angular eigenvalues for large lu is known 35] . 

In Fig. 0]and[S]we show how well the Pade function fits the data for A coming from the computer program. We also 
show the improvement compared to the second order approximation and we plot results for the real and imaginary 
parts of A. The fit parameters of the rational function representing A are shown in Table [I] for values of the angular 
parameter aa ranging between and 0.4, for the I = 2 angular mode number. 

IV. TRANSFORMATION THEORY AND THE REGGE- WHEELER EQUATION 

We have to impose some boundary conditions on the radial equation before we can obtain an equation for the 
frequencies of the quasinormal modes. We shall do this by extending the transformation theory of perturbations of 
black holes to AdS space, allowing us to construct the Regge- Wheeler- Zerilli equation for the underlying the metric 
perturbations. 

The transformation theory for Kerr black holes has been described in detail by Chandrasekhar |3l| . Following the 
Kerr case, we reparameterise the Weyl function 'Fn by 



*o = e 



„+iu>t 



(4.1) 



12 



S 

E 

CO 



10 
9 
8 
7 
6 
5 



Imaginary part of true lambda o 
Imaginary part of Pade approximation — 



9' 



10 15 
w 



20 



25 



FIG. 5: The imaginary part of the Pade' approximation of the angular eigenvalue A is shown, for aa — 0.24, / = 2 and m = 0. 



Parameter aa 


an di a2 as bi 62 


0.00 


4.0 0.0 0.0 0.0 0.0 0.0 


0.08 


4.0 0.0560 0.6068 0.0090 0.0139 0.0329 


0.16 


4.0 0.5350 0.7121 0.1039 0.1307 0.0641 


0.24 


4.0 0.9169 0.8387 0.2233 0.2177 0.1219 


0.32 


4.0 1.1610 0.9411 0.3368 0.2658 0.1770 


0.40 


4.0 1.4410 1.1728 0.3670 0.3936 0.1916 



TABLE I: This table shows the parameters we obtain by fitting the Pade approximation for A, to the numerical data coming 
out from the computer program for the calculation of the angular eigenvalues of the Teukolsky angular equation. These results 
have been obtained by fixing the angular mode number I — 2 and m = 0. 



where 



We also adopt a new r» coordinate, 



2 2 1 ~2 2,2, am X 

w = r + a = r + a + 



OJ 



ZD 

dr* — ~~^~dr 

The Teukolsky equation for Y becomes 

{ (dl + lu 2 ) + P (8* - iu) - 0} Y = 
where Q = A r u/zu 8 , and u is a polynomial 

u = (A - 6a 2 r 2 )zu 4 + 3(rA' r - A r )vu 2 - 3r 2 A r 



(4.2) 

(4.3) 

(4.4) 
(4.5) 



Part of the aim of transformation theory is to construct a function ip which satisfies the Regge- Wheeler- Zerilli 
equation 

(dl+w 2 )i, = Vi } 



(4.6) 
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with a potential function V(r). The function ip represents the metric perturbations of the black hole. 
The two radial functions Y and ip are related by an expression of the form 

^RY-^T(d,-iu)Y,= K4 (4.7) 

where R and T are functions of r and if is a constant. (In this section, we follow Chandrasekhar's notation and T does 
not denote the temperature.) The consistency of the two equations H4.4JI and (|4.6|l allows one to derive expressions 
for R, T and K, 

R = ^(u + /3A r ) (4.8) 

„ V — K 2 A r V0 2 

= u- P A r ^ 
K = k - 4u 2 f3 + 2iojK 2 (4.10) 



where 



u'A r - uA' r (4.11) 



Finding the constants k, k 2 and requires a substantial amount of algebraji^. Some help can be obtained by using 
the Starobinski constant \K\ 2 for Kerr-AdS given in |3Cf]. We obtain, 

K = \(\ + 2 + 2a 2 a 2 ) - \2a 2 a 2 (4.12) 
k 2 = ±(36M 2 + 2/3k - 12AS 2 - 10A 2 S 2 + 24Aa 2 

+24AaV - 12AS 2 a 2 a 2 - 725 2 a 2 a 2 ) 1/2 (4.13) 

/3 = ±3a 2 (4.14) 

There are four consistent solutions corresponding to the choice of signs in the above expressions. 
The hnal ingredient in the transformation theory is the potential, 

\ {K 2 TU 2 A r -v){K2W 2 U-(3v) K A r \ 

r \ tz7 4 (u + /3A r )( M -/3A,.) 2 u + [3A r P w s J { ' 1 

There are four different potentials depending on the choice of signs for k 2 and /3. Some general features of the potential 
are that it vanishes as r — > r + and approaches a non-zero constant as r — > oo (r* — > 0). The basic situation can be 
interpreted as scattering from r, ^ -oo to a boundary at r* = 0. 
The potential has poles at r = \a\ for real w < u> s , where 

amX " (4.16) 



+ a z 

(The potential can also have poles where u ± (3A r = 0.) Note that the r» variable is not a single valued function of r 
when to < u) s . This frequency range corresponds to the super-radiant regime of Kerr. For kerr-AdS, the existence of 
super-radiant frequencies can be prevented by the boundary conditions at r* = 0. 

In the axisymmetric case, m = 0, the potential is always continuous and bounded. As with Kerr, two of the 
potentials can be complex when k 2 is imaginary. 

Suppose that we are interested in the evaluation of some correlation function of a polar quantity in the conformal 
field theory which resides on the boundary of our AdS space. The AdS/CFT correspondence conjectures that this 
polar quantity can be associated to the polar metric perturbation of an AdS black hole living in the bulk. Accordingly 
we must set the axial black hole perturbation to vanish on the boundary. Since the axial perturbation function is the 
Regge- Wheeler function ip with k 2 < 0, this also means that we are led to assume that this ip goes to zero as r goes 
to infinity. 

On the other hand, we have series solutions for Y(r). Let us consider the boundary conditions on this function. 
First of all, we have from Eq. 14.71 

RY-T{d* -iw)Y->0. (4.17) 

We find, for /3 > 0, T — ► = 2ioj - (n 2 + 6Af)/(A - 6a 2 a 2 ) and R — ► a 2 A as r -> oo. The differential operator 
d rt as t — ► oo yields d Tn — ► — a 2 d x (recall x = r _1 ). If we rescale Y as 

Y = e +iur 'Z, (4.18) 
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with Z representing once again the series solution, then we find 

A^jZ + Z' = 7 (4.19) 

which represents the boundary condition for the axial problem we are considering. This equation tells us that assuming 
Dirichlet boundary condition on the Regge- Wheeler function ip, we are left with Robin conditions on the solution of 
our radial Teukolsky equation. Moreover, since 

oo oo 

Z — > ^2 a n, an d Z' — ► -r+ (n + 1) o„+i, (4.20) 

n=Q n=Q 

as x — ■* (r — > oo), the boundary conditions can be explicitly written in the form 

^ oo oo 

-^E a ™-E( n + 1 ) fl «+ 1 = ' ( 4 - 21 ) 

T + lo ° n=0 n=0 

In practice, « 2iw, and the numerical solutions described in the next section have used = 2iw, in which case 
the boundary conditions reduce to those used in |14| . 

The transformation theory which we have obtained is sufficiently powerful to give us the frequencies of some special 
quasi normal modes. The Teukolski function Y gives rise to four independent radial functions, from which we select 
ip^ with oposite signs for Ki. Working back from eq. I|4.7|) . we can relate the radial functions by 

fC^ip* = + D ± d*ip T (4.22) 

for some functions C ± and Z? ± . Now we take u> to be a root of K + (cu) = 0, and solve for %p~ to obtain a quasi normal 
mode. The frequencies are 



±(n 2 ±(nl-^n)^) (4.23) 



4/3 

Note that this is an implicit equation in the rotating case, since the constants depend on lu. In the non-rotating case, 
the frequencies reduce to uj = iK/12nM , which was found in In the rotating case, the frequencies can be pure 
imaginary or complex. These values provide a useful check on the accuracy of the numerical results presented in the 
next section. 



V. NUMERICAL RESULTS 



In this section we will present some numerical results obtained by using the procedures outlined above for the 
solution of both the radial and the angular Teukolsky equations. The radial boundary conditions used and the 
motivations for choosing them are those specified in section llVl We replace the angular eigenvalue A by the Pade 
approximation and then the quasi- normal frequencies u> are roots of eq. 14.2111 . 

Small Kerr Anti de-Sitter black holes, with the radial coordinate of the event horizon r + much smaller than the 
AdS radius oT x (r + -C oT x \ represent an unstable thermodynamic phase. Hence, we concentrate on larger black 
hole radii. In order to avoid metric singularities, the black hole rotation parameter is constrained to be less than the 
AdS radius, a < a -1 . The black hole parameters which we have used are displayed in Fig. |SJ 

The quasinormal frequencies are decomposed into real and imaginary parts, u> — lor + iu>j. Having chosen the 
metric perturbation with time dependence e +lut , we have that the pure imaginary component u>i is positive for all 
quasinormal frequencies. Fig. \7\ shows the evolution of the quasinormal frequencies of large black holes (LBH), with 
respect to the angular parameter a, for r + = 50 a -1 . Fig. [S] shows the evolution of the quasinormal frequencies of 
intermediate size black holes (IBH), with respect to the angular parameter a, for r+ = 0.8 or 1 . 

To obtain both pictures, we have used the Pade approximation for A and therefore we have been free to consider 
frequencies as large as we liked. The squares represent the results for static AdS black holes. Unfortunately, the 
numerical procedure described in section ITlI B 21 becomes unreliable for values of a higher than 0.3a -1 (for LBH) and 
0.4a -1 (for IBH), precluding us from extending the present results to very rapidly rotating holes. 

The quasinormal frequencies for rotating black holes show a behaviour which is very similar to the one found in 
the case of Schwarzschild black holes |38j | . They have a real part which displays an increase along with the imaginary 
part. The spacing between consecutive modes tends to become regular as the overtone number n gets larger, and 
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FIG. 6: The parameter ranges used for the intermediate size of black hole lie along the lines shown in the figure. The small 
circles show the parameters used in Fig. |H| 
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FIG. 7: The frequencies of the first few quasinormal modes for a Kerr black hole in Anti-de Sitter metric with r+ = 50 a . 
Results are plotted for the angular parameter aa — 0.0,0.1,0.2,0.3. The squares represent the results for Schwarzschild AdS 
black holes. The frequencies are rescaled by r+, so that Cj = lut + . 



it seems to depend on the angular parameter a. Since in AdS space the first ten frequencies are usually enough to 
ascertain their asymptotic behaviour, we have computed the quasinormal frequencies for aa — 0.1, 0.2, 0.3 in the case 
of LBH, and aa = 0.08, 0.24, 0.40 in the case of IBH, up to the tenth mode. What we have obtained for the spacing 
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FIG. 8: The frequencies of the first few quasinormal modes for a Kerr black hole in Anti-de Sitter metric with r+ = 0.8 a -1 . 
Results are plotted for the angular parameter aa — 0.0, 0.08, 0.24, 0.32, 0.40. The squares represent the results for Schwarzschild 
AdS black holes. The frequencies are rescaled by r+, so that Co = uor+. 



in the LBH regime is, 



w n+1 -uj n ~ (64.3 + 112 i)a aa = 0.1, (5.1) 
u n+ x-uj n ~ (65.2 + 113 i)a aa = 0.2, (5.2) 
uj n+ i — Lu n ~ (66.5 + 113 i)a aa = 0.3. (5-3) 



while in the IBH regime we have, 



uj n+1 - uj n ~ (1.95 + 2.11 i)a aa = 0.08 (5.4) 

u n+1 - u n ~ (1.72 + 2.46 i)a aa = 0.24 (5.5) 

u n+1 -uj n ~ (1.57 + 2.85 i)a aa = 0.40. (5.6) 

For the smaller hole, there appears to be evidence that the real part of the spacing decreases with a whilst the 
imaginary part of the spacing increases with a. From Figs. and [21 we see that the temperature is increasing with a 
in this region of parameter space. 

Both the real and imaginary parts of the frequencies become proportional to the surface gravity k for large k, which 
agrees with the behaviour for static black hole quasinormal modes ^| . This is displayed in figures and EI| m 
which we show the dependence of the real and imaginary parts of the first (n = 1) quasinormal frequency on the 
metric parameters for some values of a. Note that these pictures have been obtained by exploiting the second order 
approximation for A, and that the results for a = are those obtained by Moss and Norman in [l4| |. 

The surface gravity reaches a minimum close to the critical line in Fig. on which Cq diverges. There is some 
numerical evidence in figures |H1 and 1101 that the quasi normal frequencies are constant along the critical line, with 
u> w (3 + i)a. 



VI. DISCUSSION 



There have been investigations on AdS black holes particularly focused on their thermodynamic properties since 
the beginning of the eighties [3^ ■ Recently, Anti-de Sitter space has appeared in the surprising new context of the 
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FIG. 9: The dependence of the real part of the first non-special quasinormal mode frequency on the metric parameters is shown 
for I = 2 and a/r+ = 0.0, 0.1, 0.2, 0.3, where r+ is the radius of the event horizon, n represents the surface gravity of the black 
hole and a -1 is the Anti-de Sitter radius. 



Maldacena conjecture relating the physics of AdS in n + 1 dimensions to a strongly coupled conformal (gauge) field 
theory (CFT) in one dimension fewer. This interest has led to the study of static, charged and rotating black holes in 
AdS space in five and even seven dimensions, where the AdS /CFT correspondence can be still formulated. References 
on these investigations can be found, for instance, in [ill Ell E3 , liM E3 . EH l4(| . 

Large Kerr AdS black holes are stable against Hawking evaporation. This property is not shared by asymptotically 
flat Kerr black holes. Furthermore, Keir-AdS black holes can have rotating Killing vectors which are time-like 
everywhere outside the horizon. Following the orbits of these rotating Killing vectors, thermal radiation can rotate 
around the black hole with the same angular velocity of the hole, implying that the phenomenon known as super- 
radiance is not allowed, and the black hole is in thermal equilibrium with the thermal gas around it |42l |43| . 

Although research has been carried on for a few years in the field of quasinormal modes of black holes in AdS space 
and results have been obtained for Schwarzschild and Reissner-Nordstrom black holes, it is in this work that for the 
first time that frequencies of Kerr AdS black holes are computed. The possibility is that, through the AdS/CFT 
correspondence, such frequencies can be exploited in the calculation of the relaxation times of the perturbation of the 
corresponding thermal state in the conformal field theory on the boundary. 

It would be desirable to analyse the Kerr AdS black holes in deeper detail and see how the frequencies relate to 
the thermodynamical properties of the black holes and whether there is a direct interpretation in terms of the CFT. 
In particular, it would be interesting to look for special behaviour along the phase transition line and to examine the 
non-axisymmetric case to find out what happens to the super radiant modes. 

It is also worth mentioning that, although black hole solutions in higher dimensions than AdS 4 have been found 
|45l |46| . the calculation of their quasinormal frequencies has never been done. Therefore an extension of our work to 
AdS$ and AdS7, where the AdS/CFT correspondence can be still formulated, seems possible. 
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FIG. 10: The dependence of the imaginary part of the first non-special quasinormal mode frequency on the metric parameters 
is shown for I = 2 and a/r+ = 0.0, 0.1, 0.2, 0.3, where r+ is the radius of the event horizon, k represents the surface gravity of 
the black hole and a -1 is the Anti-de Sitter radius. The upper branch corresponds to larger values of r+ than the lower branch. 

Finally, it would be perhaps of some interest to further check if the ideas on black hole area quantisation and the 
relationship to highly damped quasinormal frequencies, can be translated into non-asymptotically flat spacetimes, as 
Anti-de Sitter space. 

APPENDIX A: RECURRENCE RELATION FOR THE SERIES SOLUTION Z 

In this appendix we describe in details how to derive the recurrence relation of the series solution Z for the radial 
Teukolsky equation associated to spin-2 perturbations of Kerr Anti-de Sitter black holes. 
Starting from equation 13 .261 for Z, we can write its derivatives with respect to x as 



& = ^a„(-xi)-"n(x-x 1 )"- 1 

71=1 
OO 

Z" = ^a n (-Xl)- n n(r ^ -l)( a ;-a;l) , 



n=2 



We re- write Eq. 13.271 as 



pZ" + (2iq - p') Z' + (6ia 2 xuj - 6Mx + 6a V - A) Z = 0. 



(Al) 



Next, consider expressions 13 . 201 and 13 . 2 1 1 for q(x) and p(x) and call C(x) and D(x) the coefficients of Z' and Z in Eq. 

ixn 



p(x)Z" + C{x)Z' + D(x)Z = 0, 



(A2) 



and let us develop all these coefficients as Taylor series around x = x\ = 1/?'+, which is the inverse of the radial 
coordinate of the event horizon of the black hole. We have, 



p{x) = a 2 + (1 + a 2 a 2 )x 2 - 2Mx 3 + a 2 x 4 , 



(A3) 



C(x) = -2iu{l + a 2 x 2 ) - p'(x) 
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(A4) 



D(x) = 6ia 2 xu - 6Mx + 6a V - A 



(A5) 



Note that p(x±) = 0. In fact p(x) = A r x 4 and x\ is, by definition, the largest root of A r (x). 

When, in Eq. IA1I we replace its coefficients with their Taylor expansions and insert into the same equation the 
series expressions for Z and its first and second derivatives, we obtain the recurrence relation, 



where the coefficients are given by 
ctn 



a n a n+1 + (3 n a n + 7„a„_i + f>„a n _ 2 = 0, 



{-xi)-( n+ V \p\xx)n{n + 1) + C{xi){n + 1)] 
1 



(A6) 



= (-si)- 



-p"{x x )n(n - 1) + C'{ Xl )(n - 1) + D( Xl ) 



\lf"{xx){n - l)(n - 2) + \c"{xi)(n - 1) + D'( Xl ) 
o 2 



In = {-X^-V 

S n = (-mr ( "- 2) [a 2 (n-2)(n-3)-4a 2 (n-2) + 6a 2 ] 
and we impose the initial conditions, a_ 2 = a_i = and ao = 1. Note that 

p'(xi) = -2(1 + a 2 xl)n 

where k is the surface gravity. 

We re-scale all the parameters with respect to x — x\, so that a = ax\ etc. The re-scaled mass 

1 



and the re-scaled surface gravity 



M = |(l + a 2 )(l + d 2 ) 



1 + 3a 2 - a 2 + a 2 a 2 



(A7) 
(A8) 

(A9) 
(A10) 



2(1 + a 2 ) 

Substituting everything into the equations for the coefficients of our recurrence relation, we finally obtain 

a n = 2(n 2 - 1)k(1 + a 2 ) + (n + l)2iw (1 + a 2 ), (All) 
/3„ = n(n - 3)(6a 2 - 6M + 1 + a 2 d 2 ) - 2(n - 3) 2iwa 2 + (6a 2 - 6M - A), (A12) 
7n = (n-2)(n-4)(2Af-4a 2 ) + (n-4)2wa 2 , 

5 n = (n-4)(n-5)a 2 . (A13) 
Notice is that, since we are considering Anti-de Sitter space with a negative cosmological constant a 2 > 0, 

(l-« 2 ) 



K > 



2(l + a 2 )' 

This is the same limit mentioned in [l4( for Schwarzschild Anti-de Sitter black holes. 



(A14) 
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